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Multiparticle State Tomography: Hidden Differences 
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We address the problem of completely characterizing multi-particle states including loss of infor- 
mation to unobserved degrees of freedom. In systems where non-classical interference plays a role, 
such as linear-optics quantum gates, such information can degrade interference in two ways, by de- 
coherence and by distinguishing the particles. Distinguishing information, often the limiting factor 
for quantum optical devices, is not correctly described by previous state-reconstruction techniques, 
which account only for decoherence. We extend these techniques and find that a single modified 
density matrix can completely describe partially-coherent, partially-distinguishable states. We use 
this observation to experimentally characterize two-photon polarization states in single-mode optical 
fiber. 
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The development of techniques for characterizing pure 
and mixed quantum states has enabled many advances in 
quantum information and related fields. Whether in or- 
der to study the effects of decoherenceU| , to optimize the 
performance of quantum logic gates to quantify the 
amount of information obtainable by various parties in 
quantum communications protocols [3[, or to adapt quan- 
tum error correction protocols to real- world situations]^], 
it is first necessary to obtain as complete a characteriza- 
tion as possible of the state of a given quantum system 
(or ensemble). In the general case of mixed states, this 
involves reconstruction of a density matrix, a mathemat- 
ically complete description of the degrees of freedom of 
interest in a quantum system. Entanglement with ex- 
perimentally inaccessible or 'hidden' degrees of freedom 
(sometimes called "the environment" ) enters the density 
matrix as a reduction in the off-diagonal coherences. In 
some quantum systems - those composed of multiple par- 
ticles that cannot be individually addressed - reduced co- 
herence can only partially describe the effects of hidden 
degrees of freedom. Another phenomenon, distinguisha- 
bility, arises when hidden degrees of freedom provide in- 
formation that could in principle be used to tell the par- 
ticles apart without necessarily leading to any changes 
in the coherences of the density matrix. This paper will 
show how a density matrix characterization of such states 
can be performed while taking into account the decoher- 
ence and distinguishability as distinct phenomena. 

Systems of particles that cannot be individually ad- 
dressed occur commonly in quantum optics and else- 
where. A central example is the Hong-Ou-Mandel 
(HOM) effect @, in which non-classical interference 
causes photon bunching at a beamsplitter. The effect 
results in photons with the same characteristics entering 
the same mode, making it impossible to individually ad- 
dress the photons, i.e., to manipulate or measure them 
individually. Many other major results in the field of 
quantum optics such as the generation of Bell states, the 
demonstration of teleportationQ, linear optics quantum 



computing [3], the generation of cluster statesQ and the 
demonstration of quantum logic gates [9l[lt3| also use non- 
classical interference and necessarily involve states with 
indistinguishable, non-individually-addrcssable photons. 

Numerous experiments have directly studied theprop- 
erties of non-individually- addressable photons [111. Il2l . 
13j . The work of Bogdanov et al. showed that the 
polarization state of two photons forms a controllable 
three-level system or qutrit suitable for many protocols in 
quantum information Q and quantum cryptography [lil 
15| and proposes a method for performing state tomog- 
raphy to characterize the density matrix of the qutrit 
state. This characterization is done under the assumption 
that the two photons in the state are indistinguishable, 
an assumption that is justified by measurement of high 
visibility HOM-likeQ two photon interference. If that 
assumption were invalid the tomography method would 
give an incorrect description of the state. 

All photons are, of course, indistinguishable in the fun- 
damental sense of obeying Bose-Einstein statistics. We 
are concerned with another, operational sense of "dis- 
tinguishable," often encountered in discussions of multi- 
photon coherence. For example, the HOM effect acts on 
photons which arrive simultaneously at a beamsplitter, 
but not on photons which arrive separated by more than 
their coherence time. We say that the photons could in 
principle be distinguished by their arrival times, and thus 
do not interfere. When we refer to distinguishability and 
indistinguishability in this paper it is this kind of distin- 
guishing information that we have in mind. Photons can 
be characterized by numerous degrees of freedom includ- 
ing arrival time, frequency, propagation direction, posi- 
tion, transverse mode and polarization. These degrees 
of freedom may or may not be experimentally accessible, 
depending on the capabilities of the experimental appa- 
ratus. In general, we describe as 'visible' those degrees 
of freedom that can be measured by a given apparatus, 
and as 'hidden' those that cannot. This paper will ex- 
plore the effect that distinguishability in hidden degrees 
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TABLE I: The measurement operators implemented in the 
tomography experiment. The detectors can detect either a 
coincidence between two photons in the H mode thereby 
implementing the projector Phh = \H1H2) {HiH2\ or a co- 
incidence between the H and V modes thereby implementing 
Phv = |.HiVa) {HiV 2 \ + \ViH 2 ) (ViH 2 \. A quarter- and half- 
waveplate at angles q and h respectively, placed before the 
detection apparatus, effectively rotate the detection operators 
to U® 2 P ({/t)® 2 where P is either Phh or P HV and U = 
exp [i% (a z cos 2h — a x sin 2ft)] exp [-£ (a z cos 2q — a x sin 2g)l , 
where a x , a y and cr z are the Pauli matrices. 

of freedom due to timing information, for instance, can 
have on measurements performed on a visible degree of 
freedom such as polarization. 

Our approach can be used on a variety of pho- 
ton states such as those created by combining non- 
orthogonal photon modes [l6| . those that combine the 
output from different spontaneous parametric downcon- 
version sources 11] and those generated by stimulated 
parametric downconversion[l7) . Although this paper 
treats systems of two photons, the same technique can 
be readily extended to a larger number of photons. This 
extension will be treated in another paper [18[. 

If one starts from the assumption that the photons 
are indistinguishable in the hidden degrees of freedom as 
in [ll[ then one will use a 3 x 3 matrix to describe a 
polarization qutrit as opposed to the 4x4 matrix used 
for distinguishable photons [191]. For a more general de- 
scription one could allow for a distinguishing degree of 
freedom and explore what limitations the fact that it is 
hidden places on measuring density matrix elements. 

A general pure state of two photons can always be 
written as a superposition of tensor products of states in 
the visible and hidden degrees of freedom 

W=E c *l^)v iB IXi> W d. (!) 

i 

where |<A;) vis , |Xi)hid are eigenstates of exchange opera- 
tors X V i S and Xhid for the visible and hidden degrees 
of freedom, respectively, with the same eigenvalue ±1. 
The requirement that the whole state be bosonic so 
that X v i s ® Xhid IV*) = IV*) guarantees that each term 
can be written with the visible and hidden parts of the 
state either both symmetric or both anti-symmetric. A 
completely general state of two photons is a mixture 
of states such as described by a density matrix 

A reduced density matrix describing only the visible 
degrees of freedom can be obtained by tracing out the 
hidden degrees of freedom in p. We define p v is = Trhid [p] 
which is sufficient to describe any measurement outcome 



on the visible degrees of freedom. Any such outcome 
has a corresponding operator B = B V i S (g> Ihid where B V i S 
acts only on the visible degrees of freedom and Ihid is 
the identity operator on the hidden degrees of freedom. 
Expectation values can be written as 

(B) = Tr [p (B vis ® Ihid)] 

= Tr vis [Tr hi d [p (B vis ® Ihid)]] 

= Tr vis [pvisBvis] , (2) 

We define projectors Ps,A onto the symmetric (S) and 
antisymmetric (A) subspaces of the Hilbert space for the 
visible degrees of freedom. Using P5 + P 4 = I and ex- 
panding p in terms of bosonic states as in eq|T] it can be 
shown that p v - ls has the property that 

p vis = P5p vis P S + PaPvisPa (3) 

This means that p v i s contains no coherences between the 
symmetric and antisymmetric subspaces, unlike distin- 
guishable particle density matrices [ljj. We now special- 
ize to the experimental situation in which the visible de- 
gree of freedom is polarization. The symmetric space is 
spanned by {\HiH 2 } , > \V1V~2)} and the antisymmet- 
ric space by \ip~) where |i/' ± ) = {\H X V 2 ) ± \ViH 2 )) /a/2. 
Pvis will be written in this basis. 

Physically the lack of coherences between symmetric 
and antisymmetric states expresses lack of information 
about the labeling of the photons. This is illustrated 
by the distinguishable-particle states I-H1V2) , IV1-H2) = 
(|V' + ) i \ip~)) /\/2, for which the phase between 
and \tp~) carries the information about the ordering of 
the photons. When this ordering is unmeasurable the 
magnitude of the coherence must be zero, as in p V j S . Nev- 
ertheless, the populations of \ip~) and are measur- 
able. For the same reason all of the coherences between 
the symmetric and anti-symmetric states in p v i s are zero. 
Thus /9 V i S divides naturally into two sub-matrices, a 3 x 3 
submatrix for the symmetric subspace and a 1 x 1 sub- 
matrix for the antisymmetric subspace as shown below: 

/ / Phh,hh Phh,4>+ Phh,vv \ \ 

_ Pip+,HH Ptp+,ip+ Pip+,VV 

Pvis — \ / 

\ Pvv,hh Pvv,i>+ Pvvyv ) 

\ (p+-,i>-) J 

(4) 

If all the population is contained in the symmetric sub- 
space then one recovers exactly the 3x3 matrix that 
would be measured under the assumption of indistin- 
guishable photons. If, however, there is population in the 
antisymmetric state then the presence of distinguishing 
information in the hidden degrees of freedom can be in- 
ferred. Distinguishing information can only be detected 
when it correlates to polarization. Measuring all the pop- 
ulation to be in \H1H2), for instance, reveals no informa- 
tion about whether the two photons are distinguishable 
because ordering and polarization are uncorrelated(22j. 
In such circumstances the distinguishing information is 
undetectable by any apparatus that is only sensitive to 
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FIG. 1: Experimental implementation of state prepara- 
tion and tomography protocols showing polarizing beamsplit- 
ters (PBS), beamsplitters (BS), non-linear /3-Barium Borate 
(BBO) crystals, a second harmonic generation crystal (SHG), 
quarter waveplates (QWP), half waveplates (HWP), polar- 
ization maintaining fiber (PMF), single-mode fiber (SMF) 
and single photon counting modules (SPCM). A spontaneous 
parametric downconversion (SPDC) crystal produces pairs of 
H and V photons. The separation between the H and V pho- 
tons is controlled with movable quartz wedges and very long 
delays can be introduced by inserting a thick piece of BBO 
into the beam. Single-mode fibre and a lOnm interference 
filter make the photons essentially indistinguishable in the 
spatio-temporal modes. Tomography is performed with a set 
of waveplates and a polarizing beamsplitter. This system can 
implement all the measurements in Table [I] 



polarization, but by the same token it is irrelevant to the 
outcome of any measurements made with that apparatus. 

By tomographic measurement of the visible den- 
sity matrix p v j s we characterize several different exper- 
imentally produced two-photon polarization states in 
polarization-maintaining single-mode fibre. As shown in 
Fig. [TJ a 50-fs pulsed TirSapph laser centered at 810 
nm was frequency-doubled to 405 nm, pumping a spon- 
taneous parametric downconversion crystal and creating 
pairs of photons. The crystal was phase- matched in a 
type-II collapsed-cone geometry[2(J so that the photons 
were not polarization-entangled, but rather emerged in 
separate H and V polarized beams. These beams were 
recombined into a single mode on a polarizing beamsplit- 
ter, creating a two-photon state. A set of matched quartz 
wedges was used to fine-tune the delay between the H and 
V photons. Also, a 3-mm-thick piece of BBO could be 
inserted anywhere in the beam path to create a large bire- 
fringent delay. This birefringent delay introduced distin- 
guishing timing information that was inaccessible to our 
detectors, which were not sensitive on the 100 fs timescale 
of the delay. This was our hidden degree of freedom. By 
controlling the delay between the pulses, the degree of 
overlap could be varied from completely overlapped to 
completely separated. Preparation waveplates inserted 
before and after the quartz wedges allowed various po- 
larization states to be created. Interference filters and 
single mode fibre served to limit distinguishing informa- 
tion in the spatio-temporal mode degrees of freedom. 

The set of projectors measured in our experiment is 
given in Table HI The detection apparatus consisted of a 
polarizing beamsplitter that projected the photons to ei- 
ther H or V and a set of single photon counting modules 
as shown in FigfT] A coincidence between detectors B 



and C measured Vhh from Table HI while a coincidence 
between A and either B or C measured Pjjv ■ These two 
fundamental measurements were rotated with the mea- 
surement quarter- waveplate and half-waveplate to imple- 
ment the full set of measurements listed in Table |TJ 

The simplest state to prepare is composed of one hor- 
izontal and one vertical photon with a variable delay be- 
tween them. When the H and V photons are overlapped 
as in FigOi, tomography shows that 98% of the popula- 
tion in the state is contained in the completely symmet- 
ric state \tp + ) indicating that the two photons are highly 
indistinguishable. When the photons are delayed by a 
time larger than their coherence time we obtain FigJ^t- 
The population splits with 45% of the population in 
and 55% of the population in \ip~), indicating that the 
photons are completely distinguishable to within the ex- 
perimental limits of our measurement. When the delay 
is less than the coherence time we obtain the state in 
FigH: with 62% of the population in and 31% of 

the population in indicating partial hidden distin- 
guishability. 

The most widely investigated state of two indis- 
tinguishable photons is the 2-NOON state [IE [Hj] 
which in terms of polarization is \2::0} HV = 

(\2 h ,0v) + \0h,2v))/V2. The state \l H ,lv) is a 
2-NOON state in the circular basis (to within a 
global phase) because of the creation operator re- 
lation a) h°> V = i (o) l,o) l — a) ro) nj /2. Experi- 
mentally, a quarter waveplate can map the circu- 
lar basis onto the H/V basis and turn |l#,ly) into 
(|2 ff ,0y) + |0H,2,/»/V2. 

When the indistinguishable state in FigJ5Ji is converted 
to a 2-NOON state in this way, all the population remains 
confined to the symmetric subspace and the generated 
state, Fig(2H, is indeed a reasonable approximation to the 
2-NOON state with a fidelity of 0.90 and a concurrence 
of 0.54. On the other hand when the same transforma- 
tion is applied to the state with hidden distinguishability 
state in Fig(2f> the state in FigJ5^ is generated. The fi- 
delity of this density matrix to the NOON state is only 
0.49, nearly identical to the 0.5 fidelity that an incoherent 
mixture of HH and VV states would have to the NOON 
state. Its concurrence is zero. In a multi-photon interfer- 
ence experiment such as that in[lj| FigJ^i would display 
high visibility interference fringes whereas the state from 
Figl2k would not. One might expect a tomography proto- 
col that assumes indistinguishable photons, such as that 
proposed in [l2^ . to break down when confronted with 
a state such as Figj2]i. To check this we used the den- 
sity matrix in Fig. [2J^ to calculate the outcomes of the 
measurements taken in [l2[ and linearly reconstructed 
an indistinguishable photon density matrix. The result- 
ing matrix (Fig. (Hf) mistakenly puts all the pop- 
ulation in the state. This density matrix will in- 
correctly predict the outcome of measurements made in 
other bases such as the diagonal basis. 

The state in Fig|2b is obtained if a 2-NOON state 
is made from indistinguishable photons and then sent 
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FIG. 2: The density matrices measured with quantum state 
tomography. The imaginary parts of (a) through (f), which 
are not shown, had all elements less than 0.05. White ele- 
ments are inaccessible to measurement, (a) indistinguishable 
H and V photons (b) distinguishable H and V photons (c) par- 
tially distinguishable H and V photons (d) indistinguishable 
photons transformed to the 2-NOON state (e) distinguishable 
photons with the same transformation applied (f) the same 
state as it would be characterized by the technique of [l| (g) 
a 2-NOON state after passage through decohering, misaligned 
BBO crystal 

through a complicated non-unitary process implemented 
by inserting a thick piece of BBO whose axis is at a small 
angle relative to the horizontal. As can be seen, the BBO 
reduces the size of the coherence between the HH and 
W terms (as well as causing some rotation), but leaves 
the \4>~) term unaffected. This demonstrates that deco- 
herence can occur without introducing distinguishability 
between the photons. Comparison with Fig. [5^ shows 



that decoherence and distinguishability are distinct ef- 
fects with different experimental signatures. 

We have developed and demonstrated experimentally 
state tomography of two-photon polarization states, in- 
cluding, for the first time, the effects of distinguishing 
information in hidden degrees of freedom. Such distin- 
guishing information destroys non-classical interference, 
is often a limiting factor in linear-optics devices such as 
quantum gates, and is not correctly described by pre- 
vious tomography schemes, which account only for de- 
coherence. Our tomography technique produces a 'vis- 
ible density matrix' which predicts the outcome of all 
polarization measurements, and which describes the ef- 
fects of both decoherence and distinguishability. This 
was demonstrated clearly with our production and mea- 
surement of a "NOON" state affected by both decoher- 
ence and distinguishability. The approach can be applied 
to other types of states and extended to larger numbers 
of particles. 
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